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Confinement-deconfinement and universal string effects 
from random percolation 

Ferdinando Gliozzi*, Marco Panero and Antonio Rago ^ 

^Dipartimento di Fisica Teorica, Universita di Torino and INFN, sezione di Torino, via P. Giuria, 1, 
1-10125 Torino, Italy. 

The 't Hooft criterion leading to confinement out of a percolating cluster of central vortices suggests defining a 
novel three-dimensional gauge theory directly on a random percolation process. Wilson loop is viewed as a counter 
of topological linking with the random clusters. Beyond the percolation threshold large Wilson loops decay with 
an area law and show the universal shape effects due to flux tube fluctuations. Wilson loop correlators define 
a non-trivial glueball spectrum. The crumbling of the percolating cluster when one periodic direction narrows 
accounts for the finite temperature deconfinement, which belongs to 2D percolation universality class. 



1. INTRODUCTION 

Understanding confinement in gauge theories 
is a major challenge of particle physics. Center 
vortices are believed to play an important role 
in explaining this phenomenon. Long time ago 
't Hooft proposed a criterion for confinement 
based on three main ingredients: assuming that 
i) there is a percolating cluster of central vortices, 
a) the Wilson loops measure the linking with the 
vortex lines and in) the vortices at different places 
are weakly correlated, the sought after area decay 
law follows. 

The flux of a vortex is conserved modulo N, 
where N is the number of elements of the center 
of the gauge group. This implies that a) the vor- 
tices are closed lines and that b) the coordination 
number of the intersection points is a multiple 
of N. While the former property is essential for 
defining topological linking, the latter does not 
take part in the argumentation. 

The unsolved difficulty in order to demonstrate 
confinement is to replace the numerical evidence 
of the weakness of the correlation among central 
vortices with a convincing proof. 

This suggests to reverse the argument using 
random percolation to define weakly correlated 
loops and regarding them as the central vortices 



of a suitable gauge theory. 

We define the following purely geometric set- 
ting: generate a sample of possible states 
{Ci,C2,...} simply by populating each of the 
links of a 3D lattice independently with occupa- 
tion probability p. The physical observables of 
this system, that we call still Wilson operators 
Wj, are associated to arbitrary loops 7 of the 
dual lattice with the following rule 

1. W.y{Ci) = 1 if no cluster of the configura- 
tion Ci is topologically linked to 7; 

2. W^{C,) = otherwise. 

The vacuum expectation value of this operator is 
defined by 



(w^) = lim yw^{a)/r 



(1) 



i=l 



* Speaker at the conference 



Note that only the closed paths of occupied links 
can contribute to W-y , while the dangling ends do 
not play any role, thus in this theory the center 
vortices have to be identified with the loops of the 
random clusters. 

If p is greater than the percolation threshold 
Pc there is a percolating cluster in the infinite 
lattice and large Wilson loops obey the area law 
by construction. 

The numerical implementation of this system 
is straightforward by comparison with usual sim- 
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ulations of ordinary gauge systems: no Markov 
process is needed to perform importance sampling 
and there are no thermalization problems and no 
critical slowing down. 

On the theoretical side, both the partition func- 
tion Z and the gauge group are trivial: Z = 1, 
G = Zi. Some obvious questions arise: \) does 
the model have a well-defined continuum limit? 
\i) do the Wilson loops exhibit the Liischer term 
and the other universal shape effects like in ordi- 
nary gauge theories? in) does the theory have a 
non-trivial glueball spectrum? iv) is it possible to 
define a fmite temperature deconfinement transi- 
tion in the same percolation picture? The answer 
to all these questions is affirmative. 

2. THE STRING TENSION 

We estimated the string tension a by fitting 
the mean values of the Wilson loops associated 
to squares of side R to the function 

{W{R)) ^ aR^ exp{-bR- a R^) . (2) 

The fits for not too small R are very good (see 
FigP). 



If the model under study has a well-behaved 
continuum limit, the scaling form of a near 
should be 

aip)^aoip~Pc)^'' , (3) 

where Pc = 0.2488126(5) on the cubic lattice 
P] and v is the correlation length critical ex- 
ponent of 3D percolation. We used the value 
ly = 0.8765(16)(2) In the range 0.258 <p< 
0.270 a one-parameter fit yields do — 8.90(3) with 
X^/d.o.f ^ 0.4. The fit is reported in Fig[21 
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Figure 1. Expectation values of square Wilson 
loops W{R) in a 64'^ lattice a.i p — 0.26 > Pc as a 
function of R. The solid line is a fit to Eq.Q. 



Figure 2. The string tension as a function of p. 
The line is a one-parameter fit to Eq.ljSJ. 



Disregarding the "string" factor R^ in {W) 
yields much worse fits. This suggests that the per- 
colation process could also account for the univer- 
sal shape corrections ascribed to effective string 
fiuctuations. 

A suitable quantity which is sensible to these 
efi'ects is W 



f = cxp(— ri^cr 



{W{L^n,L + n)) 



(4) 



{W{L,L)) 

which asymptotically (large L and L — n) should 
be, in the effective string picture, a known func- 
tion f(t) of the ratio t = ^ without any ad- 
justable parameters Fig|3|shows a nice agree- 
ment to this conjecture. 
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Figure 3. The quantity defined in Eq.l^J for three 
different values of p. The line accounts for the 
universal shape effects due to effective string vi- 
brations. No adjustable parameter is involved. 

Table 1 



Mass of the lowest glueball state 



p 


m a 




0.26 


0.2189(2) 


3.75(3) 


0.27 


0.3870(5) 


3.81(2) 



3. GLUEBALLS 

Though the occupied bond connected correla- 
tor is exactly zero by construction, the correlator 
among the occupied bonds belonging to a loop is 
non-zero, because of the constraint of being part 
of a closed path. The exponential decay rate of 
this correlator gives an estimate of the mass of the 
lowest scalar glueball state. We performed two 
numerical experiments at p = 0.26 and p = 0.27. 
The results reported in Tab^show that the devi- 
ation from the expected scaling behavior is rather 
small. 

4. DECONFINEMENT AT FINITE T 

Though the usual argument for deconfinement 
at finite T based on spontaneous breaking of cen- 
ter symmetry here is inapplicable, because there 
is no way to break Zi, it is possible to show that 



Table 2 



Percolation threshold p(r) for 1/T = 6 and 8 



1 

T 






6 


0.272355(5) 


1.497(4) 


8 


0.265615(5) 


1.510(4) 



the system under study goes through a continu- 
ous, finite temperature, deconfining transition. 

The argument runs as follows. Consider a lat- 
tice endowed with a finite temperature geometry 
L X L X 1/T with fixed p > pc and start to vary T. 
At low temperature [L ~ ^/T) there is a perco- 
lating cluster, so the system is confining. As the 
temperature increases, the system becomes more 
and more dominated by the 2D geometry. The 
crucial point is now to observe that the perco- 
lation threshold is a decreasing function of the 
space dimension. This implies that there is a 
critical value T{p) ^ 1/L where the percolating 
cluster crumbles away and the system no longer 
confines. Clearly this transition belongs to the 
universality class of 2D percolation. 

There is now an efficient Monte Carlo algo- 
rithm for studying percolation on any lattice [5]. 
It allows to calculate wrapping probabilities over 
the entire range of p in a single run. We applied 
this algorithm to evaluate the wrapping proba- 
bility around the large L directions. We chose 
L = 40,50,64,128 and l/T = 6 or l/T = 8. 
Finite size scaling allows us to extrapolate the re- 
sults to L — > oo limit, where wrapping probability 
is equal to percolation threshold. The results are 
reported in Tab[21 
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